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We extend a previously theory for the interspecific allometric scaling developed in a d + 1- 
dimensional space of metabolic states. The time, which is characteristic of all biological processes, is 
included as an extra dimension to d biological lengths. The different metabolic rates, such as basal 
(BMR) and maximum (MMR), are described by supposing that the biological lengths and time are 
related by different transport processes of energy and mass. We consider that the metabolic rates 
of animals are controlled by three main transport processes: convection, diffusion and anomalous 
diffusion. Different transport mechanisms are related to different metabolic states, with its own 
values for allometric exponents. In d = 3, we obtain that the exponent 6 of BMR is 6 = 0.71, and 
that the aerobic sustained MMR upper value of the exponent is b = 0.86 (best empirical values for 
mammals: h = 0.69(2) and h = 0.87(3)). The 3/4-law appears as an upper limit of BMR. The MMR 
scaling in different conditions, other exponents related to BMR and MMR, and the metabolism of 
unicellular organisms are also discussed. 

PACS numbers: 87.10.+e, 87.23.-n 



I. INTRODUCTION 

Several biological quantities change with organism size 
according to particular rules P, d, 13| . It is common to 
believe that these rules are related with the euclidean ge- 
ometry. However, in many cases the geometric pattern is 
not observed because physical constraints also limit how 
much an organism can be modified to cope with changes 
in scaling. Recently, considerable effort has been invested 
to understand the scaling of some of these variables un- 
der certain physical and geometrical constraints: the di- 
mensions of long bones H, the basal metabolic rate 
(BMR) [1, H,B,|g^^ and the maximum metabolic 
rate (MMR) [ll|, [U, [l| . In this paper we are nterested 
in the scaling of metabolic rate, which is the most studied 
variable in traditional allometry. 

It is accepted and empirically tested that the metabolic 
rate B and the body mass M of almost all organisms 
are connected by a power law relationship B = aM^, 
where a is a constant and b is the scaling exponent 
[11) S H, [3 • The origin and the universality of the scal- 
ing exponent of metabolic rates is a subject of great con- 
troversies and there are several debates in the literature 
[H, [H, miH, [H, [13, HH . In a recent paper [ij , we and 
a colleague proposed an unified theory for the interspe- 
cific allometric scaling of metabolism. It was developed 
in a d -f 1 dimensional space of metabolic states of organ- 
isms {d biological lengths and a physiological time). It is 
natural to include explicitly an extra temporal dimension 
in the analysis of allometric scaling because all biologi- 
cal process are time dependent. Moreover, in some cases 
this approach has produced a simple explanation for the 
problem with satisfactory results [ll|, [l^j [H, [13] ■ In that 



paper , the authors supposed that each metabolic rate 
of organisms is characterized manly by one of two trans- 
port processes, namely, convection and diffusion. In this 
paper we consider the general case in which a metabolic 
rate of 3-dimensional organisms can be characterized by 
one, two or three of the transport processes: convection, 
diffusion and anomalous diffusion. It is well known, that 
the transport in large distances is done by convection 
and the transport in small distances is done by diffusion. 
A classical example is the oxygen transported from the 
heart until the capillaries by convection and from the cap- 
illaries to the cell by diffusion. However, the mechanism 
of transport of large molecules inside a cell and between 
the cells of a tissue is still unknown and in many cases is 
suggested to be an anomalous diffusion. The three kinds 
of transport implies also that we must now deal with dif- 
ferent characteristic times. But they are all related if the 
network delivery is optimal. 



This work is organized as follows. We discuss the hy- 
potheses of da Silva, Barbosa and Silva (SBS) [ll[ and 
present a new one, and derive the main equations in Sec. 
im In Sec. mil we rederive the scenarios for BMR of 
SBS work in our present context as limiting cases. The 
BMR of mammals and birds is studied in Sec. IIVI and 
the scaling of capillaries and of aorta is obtained in Sec. 
fVl The approach to describe the MMR of endotherms is 
presented in Sec. I VII and the exponents of aorta and cap- 
illaries in the MMR conditions are discussed in Sec. IVIII 
The metabolism of unicellular organisms is discussed in 
Sec. IVIIII We summarize our results in the last section. 
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II. HYPOTHESES AND MAIN RELATIONS 

Following SBS [ll| we use the mass density 
Pd+i{Li, L2, . . . ; Ld, t) (mass per unit volume and unit 
time) and the energy density ad+i{Li, L2, ■ ■ ■ , Ld^r) 
(available energy per unit volume and unit time) to char- 
acterize the metabolic state of organisms. The use of 
energy density is justified because in the metabolic pro- 
cesses ATP cannot be supplied from outside but must be 
synthesized within the organism (within the cells). The 
efficient use of substrates by cells depends on the pres- 
ence of an adequate quantity of mitochondria as power 
house and secondly on adequate supply of fuels and 
oxygen. The fuels, which are directly related to the avail- 
able energy E, are contained inside the organism but the 
oxygen flux is supplied by outside the organism. So, en- 
ergy content is important to characterize the state of an 
organism. This can also be illustrated with the MMR 
situation. The animal runs until it has no more available 
energy. Then it fall exhausted. 

The first and second hypotheses are that natural se- 
lection enforces the constraints of scaling-invariant (inde- 
pendent of body mass) pd+i and ad+i, during evolution. 
The third hypothesis is that the scaling of the metabolic 
states is determined by the dominant dynamical trans- 
port processes of nutrients. Moreover, these processes 
are characterized by scaling-invariant quantities (diffu- 
sion coefficient, average velocity, etc.). 

Although we have d biological lengths, each one with 
its characteristic time ti (i — 1,2 d), we suppose that 
only one time r is relevant to describe the metabolic 
states. It means that the resources rates of all these 
processes must be matched ((1/ii) oc (1/^2) ■ ■ • (l/^d) « 
(l/r)) (symmorphosis principle proposed by Taylor and 
Weibel )12i]). Therefore we are considering optimal trans- 
portation networks (the new fourth hypothesis). 

It follows from the second hypothesis that E = 
ad+i rVd- Here, we have that rVd is the {d + 1)- 
volume and Vd = L1L2 ■ ■ ■ Ld- Using the power definition 
(P = dE/dt), the energy can be written in terms of the 
metabolic rate B, the power averaged over the time scale 
T, as -E = Bt. Therefore from the first and second hy- 
potheses we obtain an equation for the organism's mass, 
namely 

M = pd+i rVd , (1) 

and the following expression for the metabolic rate 

B = ad+i Vd . (2) 

Note that Eqs. ^ and ^ are valid for all metabolic 
regimes. Different metabolic scaling regimes will appear 
because there are different ways to transport nutrients. 



III. LIMIT SCENARIOS FOR THE BASAL 
METABOLIC SCALING 

A. The BMR-3 scenario 

We discuss in this section some limit cases for the basal 
metabolic rate scaling. Let us first study the BMR-3 
scenario, a lower bound for all metabolic scaling. We 
suppose that all the transportation occurs via diffusion, 
implying that 

1/2 

Li = Dit^ , with i — 1, 2 . . . , d , 

where are the scaling-invariant diffusion coefficients 
and ti are characteristic times. Since the resource supply 
rates must be matched (1/ii oc 1/^2 • ■ • 1/td oc l/r), we 
have only one time scale (r) and only one relevant length, 
namely 

Li (X L2 (X . . . (X Ld cc L . 

Note that the biological volume is given by Vd oc L'^. 
Since Li = Dir'^/'^, we obtain from Eq. ([1]) that r oc 
7\,f2/((i+2)^ This relation furnishes how L depends on M, 
namely L oc M^/('^+^\ and we can use Eq. ([2]) to obtain 
that 

L oc ^ 

T OC M^Td , 

B oc M5Td . 

In d = 3, the metabolic exponent is 6 = 3/5. Since these 
transportation processes are the slowest ones, this value 
is a lower bound for the exponent b for all metabolic situ- 
ations. Note that this scenario can, in principle, describe 
the metabolic rate of very small organisms because dif- 
fusion over short distances is fast. 



B. The BMR-2 scenario 

For larger organisms, transport by convection is uti- 
lized on large length scales because diffusion is slow. 
In the cardiovascular system of mammals, for example, 
blood circulates in a ballistic regime until the capillaries, 
where diffusion play the main role. Therefore we con- 
sider first that the BMR is driven by ballistic transport, 
namely 

L = voT , 

where the velocity vo is scaling-invariant. Then we must 
taken in account the other metabolic steps. In a "cylin- 
drical" symmetry we have Li cc L — voti (ballistic term) 

1 /2 

and d — 1 lengths Ri = Dit^ (diffusion terms), vq and 
all Di {i — 2, 3, . . . , d) are scaling-invariant. Since the 
delivery of the network is optimal {ti oc t2 ■ ■ ■ oc r), 
it follows that Ri oc DiT^^^. From Eq. (P) we ob- 
tain that T oc M2/(3+rf), implying that Li oc M2/(3+rf) 
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and Ri oc M^^^^^'^K Since the biological volume is 
Vd oc R'^~^Li, we obtain from Eq. (O that 

Li cx r oc M 3T3 ^ 

i? oc Af3T3 ^ 

l + d 

B (X M 3 + d , 

Then in d = 3, we obtain the 2/3 law. Note that this 
result was obtained without mention of the area/volume 
ratio. 



C. The BMR-1 scenario 

In this scenario all metabolic relevant lengths are re- 
lated to the ballistic transport, namely Li = Viti [i = 
1, . . . , d). Using that all characteristic times are propor- 
tional to r (fourth hypothesis), we find that Li oc r. 
In other words, there is only a single metabolic relevant 
length L oc Vqt and a single time t, both related to the 
ballistic transport. This scenario represents an upper 
bound for BMR. Since that Vd oc L*^, we find from Eq. 
©that 

L oc r oc , 
B oc . 

Therefore we find the 3/4-law for d — 3, namely B oc 
M3/4. This upper bound value is the same as these of 
West, Brown and Enquist [6| and Banavar et al. It 
is worth mentioning that Demetrius [lo| found that the b 
exponent of BMR should be in the interval [2/3, 3/4]. His 
work is based in the integration of the chemiosmotic the- 
ory of energy transduction with the methods of quantum 
statistics. 



IV. BMR OF MAMMALS AND BIRDS 

From now on we use d — 3 because cells and organisms 
are three dimensional objects. We also use a single time 
(r) for all transportation processes because all character- 
istic times are proportional to it. In order to study the 
BMR scaling of mammals and birds, let us discuss the 
nutrient transport in eukaryotic cells and between cells. 
The first biological length Li is related to the transport 
of oxygen and small molecules by diffusion, namely 

Li = Dt^'^ . 

On the other hand, large molecules can also be trapped in 
vesicles by macropinocytoses and pinocytoses and trans- 
ported in direction of the nucleus. Note that a vesicle 
can carry a relatively large quantity of fuel. Although 
the exact description of vesicular transport is still un- 
known, we suppose it as an anomalous diffusion process 
fillip, namely 

L2 = 7^,t(i/2)+- . 



The normal diffusion and ballistic transport processes oc- 
cur when X = and x = 1/2, respectively. This descrip- 
tion is supported by works |28l . |29| about the movement 
of engulfed particles on eukaryotic cells. Beads placed 
on the peripheral lamella of giant human fibroblasts are 
engulfed into the cytoplasm and move in direction of the 
nuclear region. In the lamella region the beads move 
ballistically with an average velocity ol v k, 1 fim/min 
(L oc vt). In the perinuclear region they move randomly 
within a restricted space and the authors have deter- 
mined that L oc t^^'^. Moreover, in a recent work of Neto 
and Mesquita 30] of optical microscopy, the authors con- 
clude that the movement of a macro pinosome inside a 
macrophage is ballistic {L oc vt) . The average velocity v 
varies from 0.5 ^/min to 2.0 ^m/min depending on the 
radius of the macro pinosome. Therefore is quite prob- 
able that a; w 1/2. We emphasize that a fuel vesicle is 
transported not only inside a cell but also from one cell 
to other one of the tissue. 

In the case of the BMR of mammals and birds, there 
is also a biological length 

L3 = VqT 

related to the transport by convection utilized on large 
length scales. For example, we find in mammals the car- 
diovascular system that transports blood to the capillar- 
ies. 

To obtain the exponent b we first evaluate V3 in terms 
of r, namely V3 oc DDxVqt^'^'" . Then, we use the relation 
between mass and r (Eq. [1]) to find how r depends on 
M (r oc Mi/(3+^). Finally, we use Eq. ^ to obtain how 
B depends on A/. It follows that 

B oc , 

T OC . 

The case x = give- us 6 = 2/3 and correspond to the 
BMR-2 scenario, where we have two lengths related to 
diffusion and one related to convection. This scenario 
yields the 2/3 law without mention of the area/volume 
ratio. The upper limit for the BMR, the BMR-3 scenario, 
is obtained when the three lengths are related to convec- 
tion. We obtain in this case the 3/4-law. When a: = 1/2 
the vesicular transport within a cell is ballistic and we 
have that b = 5/7 ~ 0.714. Since it is quite probable 
that X ~ 1/2, the BMR exponent of mammals and birds 
is close to 6 = 5/7. 

The empirical and predicted values of the BMR expo- 
nent b of mammals and birds are shown in Tab. HI This 
is the most analyzed and discussed allometric scaling in 
the last years. Note that all empirical values for b are in 
the predicted range, except the one (0.737(26)) obtained 
by Savage et al. [l6| with a data "binning" procedure. In 
such procedure the log-transformed data were averaged 
into equally spaced data points in order to achieve equal 
weight to all body size intervals and prevents phyloge- 
netic relatedness. However, the error bars do not exclude 
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the upper value of the range 0.714. Since their proce- 
dure has been criticized by Glazier [l^ , we note that the 
same data set without the binning procedure furnishes 
b = 0.712(13), a value in good agreement with our pre- 
diction. Perhaps the more rigorous and complete study of 
the mammalian BMR exponent is the work of White and 
Seymour [18j](see also [3l[). The authors excluded large 
herbivores of the data due to their long fast duration 
required to reach the postabsorptive state of BMR and 
obtained b ~ 0.686(14). They note that such animals are 
typically fasted for less than 72 h before the measurement 
of O2 consumption, while the postabsorptive state of ru- 
minants may require 7 days to be reached. We think that 
large mammals must be included in the data. Perhaps 
the BMR value can be obtained by a time extrapolation 
procedure, in which some measurements are realized pe- 
riodically after the initial fasting. It is quite possible that 
this procedure will slightly raise the estimation of b. In- 
teresting, the avian BMR exponent values are close to 
the lowest value of the predicted range. 

Since the heart rate scaling is obtained by using that 
F cx (X 1/t (X M~^/(^+^\ the predicted range for / 

is [—0.333,-0.286]. There is not a recent comprehensive 
analyze of this exponent for mammals. In Tab. [J we 
also shown the data for the pulse rate and respiration 
rate of mammals and birds measured in basal conditions. 
Savage et al. [16i] obtained —0.25 by binning the data of 
Brody [H (original exponent b = -0.27). Stahl [H 
claims that / = —0.25 but he not published its data. 
The empirical value for birds is more scattered and there 
is a report of an empirical value near the lower value of 
the predicted range. Although the pulse rate is more 
easy to measure than BMR, it is hard to achieve any 
conclusion about the empirical value of /. Note that 
the value —0.27 is out of the predicted range but is also 
different from —0.25. In fact, it is equidistant from —0.25 
and —0.286. It is worth mentioning that —0.27 was also 
predicted by other recent analysis of Bishop [13] ■ On the 
other hand, the respiration rate empirical exponents of 
mammals and birds have the majority of values within 
the predicted range. 



V. OTHER EXPONENTS OF BASAL 
METABOLISM FOR MAMMALS AND BIRDS 

In order to obtain other exponents, let us now char- 
acterize the network by "aorta" and "capillaries" . We 
define La, Ra and Va as the aorta length, radius and 
fluid velocity, respectively. The capillaries can be de- 
scribed by the capillary number Nc, length Ic, radius Tc 
and fluid velocity Vc- It is worth mentioning that the 
length Ic and the radius Tc of capillaries are not necessar- 
ily invariants, although, from our third assumption, we 
need some dynamical scaling-invariant quantities, like the 
blood flow speed velocity vq in the aorta or in the cap- 
illaries. The exponents related to these quantities can 
be obtained from the nutrient fluid conservation in the 



transportation network. Fluid conservation implies that 

Q ^ TTRlva = N^nrlvc , (3) 

where Q is the volume rate flow. It is clear that Q oc B is 
a natural assumption. Since Va is invariant, the aorta ra- 
dius scaling is given by Ra oc B^^"^ and the aorta length 
scaling is described by La oc VaT. These last relation- 
ships imply that the exponents qr and gl defined by 
Ra oc M'^s and La oc M°-^ are = {2 + x)/{6 + 2x) and 
ql — 1/(3 -f a;). The transition from the largest length 
scale (aorta) to the cell length scale occurs in the arte- 
rioles and capillaries. Vc must be scaling-invariant and, 
since the blood cells have the same size, we can make 
the extra assumption that rc is also scaling-invariant. 
Therefore the density of capillaries Nc/M behaves as 
Nc/M oc B/M. The capillary length can be invariant or 
mass dependent. Since the typical cell transport length 
is not scaling-invariant, the capillary length should also 
depend on M, namely Ic oc VcTc oc VcT. 

VI. THE MMR OF MAMMALS AND BIRDS 

The circulatory networks of endothermic animals make 
a transition from resting to maximum activity in such 
way that (i) the heart increases its rate and output, (ii) 
the arterial blood volume increases due to constriction 
of the veins and (iii) the total flow and muscular flow 
increase, with all muscular capillaries activated. These 
facts suggest that we have a "forced movement" during 
the characteristic time r, implying that the typical con- 
stant velocity can be written as v = agr (ao is a scaling- 
invariant acceleration). Therefore the aerobic sustained 
MMR is limited by an inertial movement accelerated dur- 
ing time T and the ballistic movement of BMR is now 
given hy L — VT — aqt^. 

In the upper limit of the MMR of animals, the MMR-1 
scenario, all lengths are related to the inertial accelerated 
movement (Li — aiT^, i = 1,2,3). Since Vis oc L"^ and 
L oc T^, we obtain from Eqs. ([T]) and ^ the metabolic 
relations: 

L (X , 
r oc Af? , 
B (X . 

This results agree with the ones obtained trough a gen- 
eralization of West, Brown and Enquist @ ideas to the 
MMR [131]. 

In the BMR description, we had i) a length (Li) related 
to diffusion of O2 and small molecules, ii) a length (L2) 
associated to the anomalous diffusion of vesicles and very 
near to a convection movement and iii) a length (L^) 
related the large scale transport of blood by convection. 
The MMR-1 scenario consider that L3 is the only relevant 
length and that Li and L2 evolved to match it. Although 
this description explains better the empirical data and is 
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more consistent with the maximal restrictions of MMR 
conditions, we must consider other cases. We call the 
MMR-2 scenario, when at least the lengths related to 
ballistic movement {L^ and L2) changes to L3 — a^r^ 
and L2 ~ a2T^, respectively. Li — Dit^^^ remains the 
same. Using Eqs. ^ and ^ we obtain that 

4 

L3 cx L2 (X M 11 , 
Li oc MTT , 
T oc M'^ , 

B (X . 

When only the length related to the large scale trans- 
ports changes (L3 — aar^, L2 — V2T, Li — £>ir^/^), we 
have the MMR-3 scenario. A similar procedure furnishes 
the following results: 

4 

L3 oc , 
L2 oc , 
Li oc , 

r oc Mi , 

B oc Mi . 

Note that the MMR and the heart frequency exponents 
should be in the intervals [7/9, 6/7] and [—2/9, —1/7], re- 
spectively. The predicted values of h and the heart rate 
exponent / for animals agree with the empirical values 
(see Tab. for animals in exercise-induced MMR con- 
ditions. But we must emphasize that the MMR data 
base is much narrower than it appears. Several refer- 
ences, for example Savage et al. [1^ and Bishop [ssj . 
represent basically the same data, namely those from the 
study of Taylor and Weibel [13] with some variation in 
the data composition. Note also that athletic species 
have a higher level of MMR than normal (non-athletic) 
ones '36']. For species with similar body mass, the MMR 
of athletic species can be 2.5 up to 5 times greater than 
the normal one. This implies that pi/di are different for 
the two groups. The MMR theory just developed must 
be valid for normal species. Since the MMR exponent 
for athletic species (6 — 0.94(2)) is very different from 
the one {h = 0.85(2)) for normal species, it is reasonable 
to assume that the inertial transport is different for the 
athletic group. If we assume that L3 — c^t^ ^ instead of 
L3 = ar^, we obtain a large exponent h = 9/10, a value 
near the empirical result. 

MMR can also be induced by exposure to low tem- 
perature. Oxygen consumption is measured during pro- 
gressive reduction of the ambient temperature until a 
decline in this consumption is observed. In these experi- 
ments, the animal loses such heat quantity that the usual 
ways to dissipate it are overwhelm. Then it is possible 
that the relevant lengths be dominate by heat diffusion 
(Li oc L2 oc L3 oc T^/^). This implies that b = 3/5. 
However, if we consider that the blood transport in the 
arterial system be also relevant we have that L3 = aoi| 



and Li oc L2 oc t^^'^. In this case we obtain that b = 3/4. 
It follows that the cold-induced MMR exponent should 
be in the interval [0.600, 0.750], in a relatively good agree- 
ment with empirical values. 



VII. OTHER EXPONENTS OF MAXIMUM 
METABOLISM FOR MAMMALS AND BIRDS 

The exponents related to the aorta are easily obtained 
from Eq. Q. Now the aorta blood velocity is not con- 
stant and grows with body mass as Va oc apT oc M^/^ 
(MMR-3). In fact this exponent should be in the range 
[0.143, 0.22]. From now on we will discuss the exponents 
always in MMR-1 scenario which is in better agreement 
with the empirical data. The exponents related to the 
aorta radius and length are given by a/j — 5/14 and 
aL = 2/7, respectively. The description of capillary scal- 
ing is not so clear. The radius Tc can be assumed invari- 
ant since the blood cells do not depend on body mass. 
There are three possibilities for the blood velocity Vc and 
the capillary length Ic'. (i) Vc and Ic are invariant, (ii) Vc is 
invariant and Ic — VcTc and (iii) Vc — cicTq and Ic — o-cT^ 
with ttc invariant. Since we are assuming that the typical 
transport length in a cell is not invariant, Ic should not 
be scaling- invariant. Let us consider the case (ii). Using 
that Tc and Vc are invariant, we obtain from Eq. ^ that 
the density of capillaries behaves as Nc/M oc M~^/^. 
If Tc oc T we obtain Ic oc M^/"^ . On the other hand, if 
Tc OC t2 and a; = 1/2 we obtain that Ic oc M^/^. 

We discuss now the aorta scaling. As already discussed 
, it is probable that a; = 1/2. In this case the BMR 
and MMR theories predict the same exponents aji and 
ul for the aorta radius and length, namely = 0.357 
and ql = 0.286. These values are in agreement with 
the empirical values (see Tab HI]) . The aorta velocity Va 
is invariant in the BMR description, in agreement with 
data (Dawson, 2003) {va oc M" °^). On the other hand, 
in the MMR description Va must depend on body mass 
{va OC Mi/'^ (MMR-1)). We do not know any empirical 
result for Va in the MMR conditions. So, it could be 
interesting to experimentally verify this simple prediction 
for normal species. 

The extra assumption that the capillaries radius is 
invariant is agreement with the empirical data 0] and 
with the theoretical-empirical estimations of Dawson (37j 
{rc oc M° °s, Ic oc M°-2i ). In the BMR description, the 
capillary velocity is also invariant. However, the capillary 
length Ic should depend on M. If a; = 1/2 we have that 
Ic oc M"-286. In the MMR description, we must also have 
Tc invariant. If we assume that Vc is invariant, we ob- 
tain that the capillary density Nc/M oc M~°-^^'^ agrees 
well with data of muscular capillary density of mammals 
(see Tab. |ll|. In this case we have that Ic — VcTc, with 
possibilities for Tc'. Tc oc t or Tc oc t2- Only the second 
possibility, together with x = 1/2 is consistent with Vc 
invariant. We obtain a result Ic oc M°-^*^ which agrees 
with the BMR description and with Dawson estimation. 
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VIII. THE BMR OF UNICELLULAR 
ORGANISMS 

We present now the BMR of unicellular organisms. 
The first biological length Li is related to the transport 
of oxygen and small molecules by diffusion, namely Li — 
Dt^/"^. The second length is related to large molecules 
that can be trapped in vesicles by macropinocytoses and 
pinocytoses and transported in direction of the nucleus. 
It is described by L2 = Dr^.T^^I'^'^^^ ^ with x = and 
X — Xjl corresponding to diffusion and ballistic trans- 
port, respectively. In order to evaluate the 6 exponent of 
unicellular organisms we need to taken in account a third 
length. A general description of this length is achieved 
by supposing that L3 = Dyt^^^"^^^^ , with y varying in 
the interval [0, 1/2] (diffusion movement: y = 0, ballistic 
one y — 1/2). To obtain the exponent 6 we first eval- 
uate ^3 in terms of r, namely Fa cx DDxDyT^^/'^'>+''+y . 
Then, we use the relation between mass and r (see Eq. 
([I]) to find how r depends on M (t oc m2/(5+2^+2!')). Fi- 
nally, we evaluate each length in terms of M and we use 
Eq. (12) to obtain how B depends on M. It follows that 

B oc M[3+2^+23,/(5+2x+2y)] ^ ^ ^j[2/(5+2x+2y)] ^ ^j^g 

case with only diffusion transport {x — Q and y = 0) give- 
us & = 3/5, a lower bound for the allometric exponent. 
The upper value for a; is a; = 1/2 (ballistic movement). 
In this case we have that h = (2 -I- j/)/(3 -f y), implying 
that the allometric exponent of unicellular organisms is 
in the interval [0.667,0.714]. 

In Tab mil it is shown the empirical values of h as 
well the predict values. The value of the exponent ob- 
tained by Hemmingsen [s^ (& = 0.75) is out of our 
predicted range. However, Prothero [33| observed that 
Hemmingsen had lumped together two metabolically dif- 
ferent groups (prokaryotes and eukaryotes). When he 
excluded bacteria, flagellates, and marine zygotes from 
Hemmingsen's data sample, he obtained h = 0.608 ±0.05 
for eukaryotes, a value just above our lower bound. On 
the other hand, Phillipson [4(3| studied the BMR scaling 
of 21 unicellular species and found h = 0.66 ± 0.092. 

IX. SUMMARY 

We developed a theory for the allometric scaling of 
metabolism based in four ad-hoc postulates: i) mass den- 
sity Pd+i and ii) available energy density Ud+i are scaling- 
invariant quantities, iii) dominant transport processes, 
which are characterized by scaling- invariant quantities, 
drive the metabolic scaling and iv) the resource rates of 
these processes are matched in order to have an optimal 
nutrient delivery. 

A lower bound for all metabolic exponents, namely 
bmin = 3/5, is found when we consider all transport pro- 
cesses as diffusive ones. 

The BMR of mammals and birds is obtained when we 
have A) diffusion, describing the transport of oxygen and 
small molecules, B) ballistic transport (L3 oc vqt)^ which 



is related to the blood delivery in large scale and C) 
anomalous diffusion that represents the vesicular trans- 
port inside a cell and between cells of a tissue. Assuming 
that the last process is very close to the ballistic transport 
we obtained that 6 = 5/7. This value is in good agree- 
ment with the best empirical estimation for BMR (0.69), 
obtained by White and Seymour (2005) for mammals 
without ruminants. We believe that these large mammals 
should be included in some way in the empirical analy- 
sis, implying that the real empirical h value for mammals 
should be around 0.71. The 2/3-law is obtained when 
the anomalous diffusion process is near normal diffusion. 
Therefore 6 = 2/3 is a lower bound for BMR of mammals 
and birds. On the other hand, the 3/4-law appears as an 
upper bound for BMR since this value is obtained when 
all transport processes are ballistic. Therefore the h ex- 
ponent for mammals and birds is the interval [2/3, 3/4]. 
Interesting, the empirical value for birds is close to the 
low limit. The predict interval fort the exponent related 
to heart rate (or respiration rate) is [—1/3,-1/4] and the 
most probable value is —2/7. They are in agreement with 
the empirical values for mammals and birds. However, 
this exponent was not recently studied in a comprehen- 
sive way as was done with the h exponent. 

The aerobic sustained MMR is described by an iner- 
tial movement accelerated during time r (L3 oc apT^). 
The upper limit for the MMR exponent {h = 6/7) is ob- 
tained when all transport processes are proportional to 
the accelerated one and the lower limit {h = 7/9) occurs 
when only the length related to the large scale transport 
(£3) changes to the accelerated movement. Since during 
strenuous exercise the transport systems are stressed to 
their uttermost, we believe that the upper limit describe 
better the MMR scaling. The predict value for 6 and for 
the heart rate exponent are in good agreement with the 
empirical values. However, the data base is still narrow. 
The cold-induced MMR is studied by considering that 
usual heat transport processes are overwhelm and that 
we have a new metabolic state where only heat diffusion 
is important. The different empirical exercise induced 
MMR exponents obtained for athletic species and non- 
athletic one can be explained qualitatively by assuming 
that the accelerated movement for athletic species is dif- 
ferent (L3 oc cqt^ , for example). 

The exponents related to the aorta and capillaries of 
mammals are obtained through fluid conservation. Aorta 
blood velocity is scaling-invariant in BMR conditions but 
grows with mass in the exercise-induced MMR situation. 
This exponent, which is predicted to be in the interval 
[0.143, 0.22], was never measured. The empirical deter- 
mination of this exponent seems to be easy and interest- 
ing. Moreover, it can be a experimental test of the impor- 
tance of the transportation processes for the metabolic 
scaling. The exponents characterizing the length and 
radius of aorta and capillaries have the same values in 
the BMR description and in the upper limit of exercise- 
induced MMR situation. The predict values agree with 
the empirical ones. On the other hand, the capillary 
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density must be described by the MMR scenario and the 
predict value also agrees with the experimental value. 

Finally we discussed the BMR of unicellular organisms. 
In this case two transport processes are diffusion, related 
to oxygen and small molecules transport, and anomalous 
diffusion, related to vesicular transport. The third trans- 
port process was described also as an anomalous diffusion 
one because it allow us to change the transport mecha- 
nism from diffusion until a ballistic movement. The pre- 
dicted range [2/3, 5/7] was compared with the few results 
of the literature. 
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TABLE I: AUometric exponents related to the basal 
metabolism for mammals and birds. Under parenthesis is 
the error in the last significative of the observed quantities. 



Predicted 


Observed 


Comment 


Ref. 


Basal metabolic rate 


0.714 


0.712(13) 


mammals 


fl6] 


[0.667,0.714] 


0.737(26) 


binning data 


fl6] 


0.666 (BMR-2) 


0.668(25) 


small mammals 


fl5] 


0.750 (BMR-1) 


0.710(21) 


mammals 


im 




0.664(15) 


birds 


\M 




0.686(14) 


large mammals excluded 


fl8] 




0.669 


birds 


[41] 


Heart rate 


-0.286 


-0.25(2) 


mammals, unknown data 


[33] 


[-0.333, -0.286] 


-0.25 


mammals 


mj 


-0.333 (BMR-2) 


-0.27 


mammals 


[32, 43, 44] 


-0.250 (BMR-1) 


-0.25(3) 


mammals, binning data of 


[16J 




-0.26 


mammals 


[45] 




-0.27 


mammals 


[46] 




-0.23 


birds 


[47] 




-0.209 


birds 


[48] 




-0.33(6) 


birds 


[49] 


Respiration rate 


-0.286 


-0.260(5) 


mammals 


[33] 


[-0.333, -0.286] 


-0.28 


mammals 


m 


-0.333 (BMR-2) 


-0.28 


mammals 


m 


-0.250 (BMR-1) 


-0.25 


mammals 


[51J 




-0.26(6) 


mammals, binning data of 


[16] 




-0.31 


birds 






-0.20 


birds 


[48J 




-0.33 


passerines 


[42] 




-0.28 


nonpasserines 


[42] 
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TABLE II: AUometric exponents related to maximum 
metabolism and other exponents for mammals and birds. Un- 
der parenthesis is the error in the last significative of the ob- 
served quantities. 



Exercise-induced maximum metabolic rate 


0.857 (MMR-1) 


0.828(70) 


mammals, binning data 


[16] 


[0.778, 0.857] 


0.88(2) 


standard animals 


[35] 




0.872(29) 


mammals 


[36] 




0.94(2) 


athletic 


[36] 




0.85(2) 


non- athletic 


[36] 




0.882(24) 


marsupials 


[53] 




0.87(5) 


mammals 


\M 


Heart rate 


-0.143 (MMR-1) 


-0.16(2) 


mammals 


[34] 


[-0.222,-0.143] 


-0.15 


mammals 


[45] 




-0.17(2) 


birds, flight 


[49] 




-0.146 


birds, flight 


[48] 




-0.157 


birds 


[48] 


Cold-induced maximum metabolic rate 


[0.600, 0.750] 


0.65(5) 


mammals 


m 




0.772(30) 


marsupials 


m 




0.789(40) 


theria 


[53] 


Capillary density 


-0.143 (MMR-1) 


-0.14(7) 




[2, 54] 


Aorta radius 


0.357 (BMR and MMR-1) 


0.36 


mammals 


[2,46] 


[0.333,0.357] (BMR) 


0.335 


mammals 


[44] 


Aorta length 


0.286 (BMR and MMR-1) 


0.32 


mammals 


[2, 44] 


[0.286,0.333] (BMR) 


0.31 


mammals 


[46] 



TABLE III: BMR allometric exponent 6 (B ~ M'') for uni- 
cellular organisms. Under parenthesis is the error in the last 
significative of the observed quantities. 



Exponent b 


Predicted 


Observed 


Ref. 


0.714 


0.75 


[38] 


[0.667, 0.714] 


0.66(9) 


[40] 


0.600 (MR-3) 


0.608(5) 


[39] 



